Electronic structure of heavily-doped graphene: the role of foreign atom states 
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Using density functional theory calculations we investigate the electronic structure of graphene 
doped by deposition of foreign atoms. We demonstrate that, as the charge transfer to the graphene 
layer increases, the band structure of the pristine graphene sheet is substantially affected. This is 
particularly relevant when Ca atoms are deposed on graphene at CaCe stoichiometry. Similarly to 
what happens in superconducting graphite intercalated compounds, a Ca bands occurs at the Fermi 
level. Its hybridization with the C states generates a strong non-linearity in one of the vr* bands 
below the Fermi level, at energies comparable to the graphene E2g phonon frequency. This strong 
non-linearity, and not manybody effects as previously proposed, explains the large and anisotropic 
values of the apparent electron-phonon coupling measured in angular resolved photoemission. 

PACS numbers: 71.15.Mb, 74.25.Jb, 79.60.-i 



The discovery of superconductivity in CaCeiE Hi Hi 
demonstrates that superconducting critical temperatures 
(TcS) as large as 15 K can be obtained in Graphite In- 
tercalated Compounds (GICs). Recently, deposition of 
foreign atoms onto a graphene monolayer [J| has been 
achieved. This finding could considerably widen the 
number of GICs, since the constraints for deposition onto 
graphene are milder than those required for graphite in- 
tercalation. Consequently even larger TcS could be dis- 
covered in these systems. 

So far, the electronic structure of atoms deposed on 
graphene has been interpreted in terms of the pristine 
graphene band structure 0, @, @1- This is correct 
if the charge transfer to the graphene layer (doping) is 
weak. However, this is questionable for larger dopings, 
since foreign atom states could affect the band structure 
below the Fermi level (e/)- Indeed in GICs, an inter- 
calant band crosses e/ 0, [3, 11 1 . As far as K deposi- 



tion on graphene is concerned, detailed Angular Resolved 
Photo Emission ( ARPES) measurements p, 0, 0, Q have 
shown that (i) the graphene band-structure below e/ is 
weakly affected by the presence of K atoms, (ii) a marked 
kink occurs at 0.195 eV below ej, an energy correspond- 
ing to the E2g graphene phonon frequency and (iii) the 
electron-phonon coupling extracted from ARPES mea- 
surements is 5.5 times larger than what could be ex- 
pected on the basis of a rigid doping of the graphene 
bands [l^, 12 1. In a very recent paper [l^l deposition 
of Ca atoms on graphene at doping as large as that of 
CaCg was achieved and, most surprising, a massive en- 
hancement of the electron-phonon coupling was reported. 
The apparent electron-phonon coupling has been shown 
to be very anisotropic with values ranging from 0.5 to 
2.3, suggesting that an heavily-doped graphene mono- 
layer could be superconducting with large critical tem- 
peratures and large anisotropic superconducting gap(s). 
Such large and anisotropic values of the electron-phonon 
coupling were interpreted as due to many-body effects 



and to the occurrence of a Van hove singularity. However 
neither measurements nor calculations of the electronic 
density of states were performed. 

In this work we use density functional theory cal- 
culations to interpret the huge and highly anisotropic 
electron-phonon coupling observed in heavily doped 
graphene by ARPES in terms of the CaCe monolayer 
band structure. 

In ARPES the spectral weight is measured, namely, 
for a given band index. 
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where I]g^jj(k, e) and I]g^jj(k, e) are the real and imaginary 
parts of the electron self-energy I]an(k, e), and ek are the 
single particle bands. The total electron self-energy in- 
cludes contributions from all the interactions in the sys- 
tem. The electron-phonon coupling parameter is: 
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where I](k, e) is the electron-phonon contribution to the 
electron self-energy and S (k, e) its real part. In what 
follows we assume that the dominant contribution to the 
electron self-energy is given by the electron-phonon cou- 
pling, Eaii(k, e) « S(k, e). 

The maximum position in the spectral weight is given 
by the relation e^'^^ — Ck ~ 5] (k, e'^'^^) = 0. Linearizing 
S (k, e) « — AkC leads to Ak = ^fer — 1- If the bare bands 
behave linearly , namely ek = hv^k, and if th e renormal- 
ized can be linearized at ef, e^^^ = hvfk, [12, llJl then Ak 
and 

Ak = — 

where, w/ is obtained from a linear fit to the maximum 
position in momentum distribution curves (MDCs) at en- 
ergies close to e/, while vq is obtained from a linear fit in 
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an appropriate energy window below the kink where the 
bare bands are hnear. Note that if the bare bands are 
not hnear then Eq. [3] is incorrect. 

In experiments a finite resolution affects substantially 
the energy window close to e/ so that the determination 
oi Vf is non-trivial 
then computed: 



15[. The the following quantity is 
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where now uarpes is obtained from a linear fit to the 
maximum position in MDCs in an energy window rang- 
ing from the kink energy up to the Fermi level. Clearly 
Aarpes ~ A only if Vf w warpes- Consequently the en- 
ergy window chosen to fit warpes is crucial. This has 
been demonstrated to be important in graphene at low 
electron-doping [l2|. Indeed, in the case of rigid band 
doping of the graphene tt* bands, where the electron- 
phonon coupling can be calculated analytically 13, 1^ . 
it has been found that Aarpes ~ 2.5A. Thus, in this case, 
Aarpes cannot provide a quantitative information of the 
electron-phonon coupling since warpes depends critically 
on the fit procedure. 

When large dopings are considered, as in the case of 
a CaCe monolayer in ref. [l^l, the commonly accepted 
interpretation of the band structure in term of rigid dop- 
ing of graphene tt* states becomes questionable. For this 
reason we calculate the single particle band-structure of 
monolayer and bulk CaCe using density functional the- 
ory (DFT). We assume that the Ca deposition leads to 
an ordered x \/3 structure with axes rotated of 30° 
respect to the standard C2 crystal structure. Simulating 
a single layer requires very large cells to converge due to 
the finite electric dipole formed by the Ca donor and the 
graphene acceptor. The problem can be solved using a 
double slab geometry with two specular layers of CaCg 
separated by vacuum along the z direction. This geom- 
etry has zero net electric dipole. The interstitial space 
between Ca atoms on different layers is 11 A while that 
between the two Graphene layers is 10 A. We optimize 
the in-plane graphene lattice parameter (a) and the dis- 
tance (z) between the Ca atoms and the nearby graphene 
layer, obtaining a = 2.47 A and z = 2.315 A. Note that 
in bulk CaCg structural optimization leads to a = 2.50 
A and z — 2.60 A 11 1. Electronic structure calcula- 
tions are performed using the espresso codeflG^ and the 
generalized gradient approximation . We expand the 
wavefunctions and the charge density using a 35 Ry and 
a 600 Ry cutoffs, respectively. The electronic integration 
for the double slab geometry has been performed using a 
8x8x8 k-points mesh and a Gaussian smearing of 0.05 
Ryd. For graphene we use a 30 x 30 x 2 k-points mesh. 

In fig. [T]we show the band structure of CaCg mono- 
layer compared to that of graphene using the same in- 
plane lattice parameter, a = 2.47 A. The energies of the 
Dirac point in the two structures have been aligned. In 
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FIG. 1: (Color online) (a) Band structures of CaCe Mono- 
layer (continuous line) and graphene (red-dashed) using the 
same in-plane lattice parameter, namely a = 2.47A (dashed). 
The size of the circles indicates the percentage of Ca char- 
acter (see ref. [l3|) in a given band, (b) Fermi surface of 
CaCs monolayer. The special points labels refer to the CaCe 
hexagonal Brillouin zone. 



the CaCg monolayer Brillouin zone the Dirac point is 
refolded at F due to the -x/S x v^RSO" surface deposi- 
tion. Because of this superstructure, a small gap opens 
at the Dirac point in CaCg monolayer. The Fermi level 
is crossed by two Carbon tt* bands and a third band 
mainly Ca in character. The steepest tt* band (labeled 
TT^*) is very similar to that of C2, while the other tt* (la- 
beled TTj) is substantially affected by the Ca recovering. 
The most striking feature is, however, the occurrence of 
an additional band at ej, having dominant Ca charac- 
ter. This band is the so called "intercalant" band in 
GICs [3, 113] , descending to lower energies respect to its 
position in CaCg bulk. The bottom of the intercalant 
band is at —1.15 eV respect to e/ in bulk CaCg and 



3 




FIG. 2: (color online) (a) Band structure of CaCs Monolayer (continuous and dotted lines) and of Graphene (dashed lines) 
plotted in the Graphene Brillouin zone, (b) Fermi surface of CaCe as a function of doping. Square and diamonds are data 
from ref. and refers to the dopings labeled Ca (top) and Ca (under) K 2 x 2 (top). Circles are theoretical calculations. 
The size of the circles is proportional to Aarpes- (c) Polar plot of Aarpbs- The high-symmetry points are labeled respect to 
the G2 Brillouin zone. 



is now at —1.65 eV in the monolayer. Correspondingly 
the graphene Dirac point is at « —2.04 eV in the bulk 
and at « 1.39 in the monolayer, in nice agreement with 
experiments 13[ . This effect is well known; in alkaline- 
earths GICs, the intercalant band descends as the inter- 
layer distance is reduced [111, IHl • 

Finally, in CaCg monolayer there is a marked avoided 
crossing between the ttJ band and the Ca-band along 
TM. In this direction, approaching the M-point, there 
is a large hybridization between the two bands, as it is 
evident from the change in the Ca character of the band. 
Most important, a large deviation from linearity occurs 
in the ttj band at an energy of « 0.2 eV below the Fermi 
energy. This energy is very close to the E2g phonon fre- 
quency [l9| of undoped graphene at the T point, namely 
UJE,,T = 0.195 eV. 

The Fermi surface of CaCg monolayer is shown in fig. [T] 
(b) in the CaCg monolayer Brillouin zone and in fig. [2] (b) 
in the C2 Brillouin zone. In fig. [T] (b) the Fermi surface 
is composed of a 6-points star centered at F formed by 
the intersection of two triangular shapes. In addition 
there are 6 hole pockets centered around K an K'. The 
hexagon within the star is due to the ttI band, while 
the outer perimeter of the star is composed by the ttJ 
band hybridized with Ca states. Along TK the star has 
dominant C character, while along TM it has strongly 
hybridized Ca and C character. The hole pockets are 
mostly Ca. In fig. [2] (b) one of the triangular shapes 
composing the star is compared with the Fermi surface 
determined by ARPES [13|. A very good agreement is 
found for the surface at the highest experimental doping. 

In fig. [21 the CaCg monolayer and C2 bands are plotted 
in the C2 Brillouin zone. The aforementioned deviation 
from linearity in the ttj bands, now occurring along MK, 
results in an apparent kink. However this kink is not 



due to the electron-phonon interaction or to many body 
effects, but it is due to the non-linearity of the single 
particle bands. The change in slope is induced indirectly 
by the hybridization with the intercalant band. 

In order to compare with experimental data, we esti- 
mate the "apparent" electron-phonon coupling Aarpes 
(Eq. [4|) generated by the change in slope of the single 
particle band. We obtain vq and uarpes in Eq. [3] by 
linear fits to the band at several directions departing 
from the K point in the C2 Brillouin zone. The fits to 
obtain warpes and vq are performed in the energy range 
-0.2eV < e - 6/ < OeV and -l.OeV < e - e/ < -0.2eV, 
respectively. A strong anisotropy in Aarpes is obtained 
as a result of the single-particle band-structure. Values 
as large as 2.5 are found along the KM direction. Along 
TK, where no kink is present in the bare bands, a zero 
value is found. 

To compare with experiments the contribution due 
to the electron-phonon coupling in the CaCe monolayer 
should be added. In ARPES spectra in ref. [iSj the only 
visible contribution is that due to the in-plane graphene 
phonons at ~ 0.2 eV. In bulk CaCe this contribution is 
0.11 [13] ■ As we have shown in ref. [l3|, the fitting pro- 
cedure used in ref. [l3| overestimates A of a factor 2.5. 
Thus to the anisotropic bare band contribution we added 
an isotropic electron-phonon contribution of 0.275. 

The results are illustrated as polar plots in Fig. [2] (b) 
and (c) and are compared to the experimental data of ref. 
[ist in Fig. [21 (c). We reproduce the strong anisotropy 
observed in ARPES experiment [l3| and the huge val- 
ues of Aarpes along the KM direction. In the KF di- 
rection, where the bare band contribution is zero, the 
value of Aarpes is given entirely by the electron-phonon 
contribution. In this direction the values of A for bulk 
CaCg underestimates the experimental one. This could 
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be due a larger density of states at the Fermi level in 
CaCe monolayer respect to bulk CaCe- Indeed we found 
that the monolayer density of states at the Fermi level is 
1.7 time larger than that in the bulk. 

In this work we have studied the electronic structure 
of foreign atoms deposed on graphene. We have shown 
that as the charge transfer to graphene becomes signifi- 
cant, the band structure cannot anymore be interpreted 
in terms of the pristine graphene it* bands. When Ca is 
deposed on graphene at stoichiometry comparable with 
CaCe, a Ca band occurs at the Fermi level. This band 
is the surface analogue of the intercalant band in su- 
perconducting graphite intercalated compounds 0, [lo| , 
that is crucial to understand the superconducting be- 
havior of these compounds. The Ca bands strongly hy- 
bridize with the C tt* states and induces a marked non- 
linearity in one of the tt* bands at energies compara- 
ble to ef — CL;E2gr- This non-linearity explains the large 
and strongly anisotropic values of Aarpes in CaCg mono- 
layer. However, as shown in this work, these large and 
anisotropic values do not provide any information on the 
real electron-phonon coupling in the system. 
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